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Results on Dalitz analysis of three-body charmless B + — > K + -k + -k~ and B + — * K + K + K~ 
decays are reported. We also present preliminary results on the studies of direct CP violation 
in three-body decay B — > K ± tt + -k~ . The analysis is performed using a large data sample 
collected with the Belle detector. 



1 Introduction 



Studies of B meson decays to three-body charmless hadronic final states are a natural extension 
of studies of decays to two-body charmless final states. Multiple resonances occurring nearby in 
phase space will interfere and a full amplitude analysis is required to extract correct branching 
fractions for the intermediate quasi-two-body states. In addition to being a rich laboratory for 
studying B meson decay dynamics, three-body charmless hadronic final states may provide new 
possibilities for CP violation searches^. 



2 Apparatus, Data Sample & Event Selection 

This analysis is based on a 140 fb _1 data sample (152 million BB pairs) used for Dalitz analysis 
and on a 253 fb _1 data sample (274 million BB pairs) used for searches for direct CP violation in 
decay B^ — > K^tt+tt^ . The data are collected with the Belle detector operating at the KEKB 
asymmetric-energy e + e~ collider with a center-of-mass energy at the Y(4«S) resonance. 

The Belle detector^! is a large-solid-angle magnetic spectrometer based on a 1.5 T super- 
conducting solenoid magnet. Charged particle tracking is provided by a three-layer silicon 
vertex detector and a 50-layer central drift chamber (CDC) that surround the interaction point. 
Charged hadron identification is provided by dE/dx measurements in the CDC, an array of 




1188 aerogel Cerenkov counters (ACC), and a barrel-like array of 128 time-of-flight scintillation 
counters (TOF); information from the three subdetectors is combined to form a single likelihood 
ratio, which is then used for pion, kaon and proton discrimination. For charged kaon identifi- 
cation, we impose a requirement on the particle identification variable that has 86% efficiency 
and a 7% fake rate from misidentified pions. Charged tracks that are positively identified as 
electrons or protons are excluded from the analysis. 

The dominant background to studied processes is due to e + e~ — > qq (q = u, d, s and c 
quarks) continuum events. This type of background is suppressed using variables that charac- 
terize the event topology. This allows to reject about 98% (92%) of the continuum background 
in the B ± -> K^+tt' (B ± -> K^K+K') decay while retaining 36% (70%) of the signal. A 
detailed description of the continuum suppression technique can be found in Ref.^ and refer- 
ences therein. The background that originates from other B meson decays is studied using a 
large sample of Monte Carlo generated events. We find that the dominant BB background in 
the K+TT+TT- final state is due to B + -> D°tt + , D° -» A+vr" and also B + -> J/i/;(i/)(2g))K + , 
J/t/j(tp(2S)) — > [i + fJ>~ decays. We veto these events by requiring \M(Kn) — Mn\ > 0.10 GeV/c 2 , 
|M(vr + 7r-) M+At - - Mj/^l > 0.07 GeV/c 2 and |M(7r + 7r-) M+/1 - - M^ 2S) \ > 0.05 GeV/c 2 , where 
subscript indicates that the muon mass assignment was used for charged tracks to calculate 

the two-particle invariant mass. We also remove B + — > D°K + , D° — ► it + tt~ signal by requiring 
\M(7r + ir~) — Mjj\ > 15 MeV/c 2 . To suppress the background due to tt/K misidentification, we 
also exclude candidates if the invariant mass of any pair of oppositely charged tracks from the 
B candidate is consistent with the D° — ► K + tt~ hypothesis within 15 MeV/c 2 , regardless of the 
particle identification information. The most significant background from charmless B decays 
is found to originate from B + — > rj'K + followed by rj' — > 7r + 7r _ 7, B + —* p°ir + and B° — > K + ir~ 
decays. These backgrounds cannot be removed and are taken into account when fitting the data. 

The dominant background to the K + K + K~ final state from other B decays is due to 
B — > Dh process, where h stands for a charged pion or kaon. To suppress this background, 
we reject events where any two-particle invariant mass is consistent with D° — ► K + K~ or 
D° — > K + n~ within 15 MeV/c 2 regardless of the particle identification information. We find no 
charmless B decay modes that produce a significant background to the K + K + K~ final state. 

3 Results on Dalitz Analysis of B + -> AT+vr+vr and B+ -> K + K + K~ 

Results of this analysis are described in detail in Ref. ^. There are 2584 A ± 7r + 7r _ and 1400 
K^K + K~ events in the B signal region that satisfy all the selection requirements. Two-body 
invariant mass distributions for these events are shown in Fig. ^ We find that the best fit to 
events is obtained with a matrix element that is a coherent sum of the A* (892) ir + , 
A^(1430)°7r+, p(770)°A+, / (980)A+, / X (1300)A+, XcoK + quasi-two-body channels and a 



Table 1: Summary of branching fraction results. The first quoted error is statistical, the second is systematic and 

the third is the model error. 



Mode 



K ± tt + tt charmless total 
K*(892)°7r+, K*(892)° -> K+ir- 
K *(1430) vr+, #5(1430)° -> K+tt- 
K*(1410)°7T + , K*(1410)° -» K+tt- 
K*(1680)°vr+, K*(1680)° -» K+ir- 
K 2 *(1430) tt+, K 2 *(1430)° -» IT+tT 
p(770)°i^ + , /?(770)° vr+TT- 
70(980)^+, / (980) ^+^- 
/ 2 (1270)K+, / 2 (1270) - TT+vr- 
Non-resonant 



£(5+ - 



Rh+) 
hr) x 10 6 



,6(5+ -» Rh + ) x 10 b 



46.6 ±2.1 ±4.3 



+0.38 



-0.57 



6.55 ±0.60 ±0.60 
27.9 ± 1.8 ± 2.6±|;| 

< 2.0 

< 3.1 

< 2.3 

4.78 ± 0.75 ± 0.44+°;^ 
7.55 ±1.24 ±0.69l^ 

< 1.3 



9.83 ±0.90 ±0.90 



+0.57 
-0.86 



45.0 ± 2.9 ± 6.2± 1 |;J 



4.78 ± 0.75 ± 0.44±g;|^ 



17.3 ± 1.7 ± 1.6 



K ± K + K charmless total 
(f)K + , (f> -» K+K~ 
0(168O)K + , 0(1680) -» K+if- 
/ (980)^+, / (980) -» K + K~ 
&(1525)K+, $(1525) - 
a 2 (l320)K + , a 2 (1320) -» iT+iT 
Non-resonant 



+17.1 
-7.8 



4.72 ± 0.45 ± 0.35iaia 

< 0.8 

< 2.9 

< 4.9 

< 1.1 



30.6 ±1.2 ±2.3 
9.60 ± 0.92 ± Q.ntoM 



1.37 ±0.28 ±0.121^ 
0.86 ± 0.26 ± 0.06±g;g§ 



24.0 ± 1.5 ± ljjrjjj 



XcoK + , XcO -> 7T+7T 
Xc0^ + , XcO - 
Xc0-^ + combined 



196 ± 35 ± 33l^ 7 



non-resonant component. A channel fx(^300)K + (with mass and width of .fx (1300) to be 
determined from the fit) is added to account for an excess of signal events visible in M(ir + Tr~) 
spectrum near 1.3 GeV/c 2 . Results of the best fit are shown in Figs. ^ (a,b). The mass and 
width of the /x(1300) state obtained from the fit are consistent with those for the /o(1370), 
however more data are required for more definite conclusion. To test for the contribution of other 
possible quasi-two-body intermediate states such as J PC*(1410) vr + , K* (1680)° vr+ , iT 2 *(1430)°7r + 
or / 2 (1270)iT + , we include an additional amplitude of each of these channels in the matrix 
element one by one and repeat the fit to data. None of these channels have a statistically 
significant signal. Branching fraction and upper limit results are summarized in Tabled 

The best fit to K :iz K + K~ events is obtained with a matrix element that is a coherent sum of 
the (j)K + , /x(1500)i<T + , XcoK + quasi-two-body channels and a non-resonant component, where 
the /x(1500)-?T + (with mass and width of /x(1300) to be determined from the fit) channel is 
added to describe the excess of signal events visible in K + K~ mass spectrum near 1.5 GeV/c 2 . 
As there are two identical kaons in the final state, the decay amplitude is symmetrized with 
respect to interchange of two kaons of the same charge. Results of the fit are shown in Fig. ^ (c) 
and summarized in Table In order to check the sensitivity of the data to the spin of the 
/x(1500) state, we replace the scalar amplitude by a vector or a tensor amplitude for the 
/x(1500): the scalar hypothesis gives the best fit. To estimate the possible contribution from 
other quasi-two-body channels, we include an additional decay channel and repeat the fit to the 
data. In particular we test the c/>(1680)X + , / 2 (1525)-KT + and a 2 (1320)i , T + channels. In all cases 
the fit finds no statistically significant signal for the newly added channel. Since we observe a 
clear /o(980)/C + signal in the K + tt + it~ final state, we try to include the fo(980)K + amplitude 
in the B + —* K + K + K~ matrix element as well: no statistically significant contribution from 
this channel is found. 




Figure 2: 



M(tt + tt ) mass spectra for B (left) and B + (right). Points with error bars are data, the open 
histogram is the fit result and the hatched histogram is the background component. 



4 Search for Direct CP Violation in — > 



TT 1 TT 



In the analysis of the direct CP violation in — > K ± tt + tt~ decays we use larger data sample 
of 253 fb^ 1 . To search for direct CP violation we make a fit to K~tt~tt + and K + tt + tt~ samples 
separately. In the fit to B~ and B + subsamples we fix all the parameters not sensitive to CP 
violation (such as masses and widths of resonance states) at values determined from the fit to 
the combined B sample. Results of calculations of CP violation asymmetries defined as 



Act 



B{B- 



(1) 



are given in Table |2J In all except for B^ — > p(770)° K channels the measured asymmetry 
agrees with zero within la — 1.5a (after the model uncertainty is accounted). A 2Aa hint for 
a large direct CP violation is observed in the decay B^ — > p(770)° . This is demonstrated 
in Fig. [21 where the 7r + 7r~ invariant mass distribution for the p(770) — /o(980) mass region is 
shown. Analysis with a large data sample is required in order to check this result. 



Table 2: Results of .Acp calculation. Quoted errors are statistical, systematic and the model related error. 



Mode 



K ± ir + ir Charmless 
K*(892)°vr + , K*(892)° -► K + tt- 
K^(1430)^ + , JfJ(1430)° -> K+7T- 
p(770)°K + , p(770)° -» iT+iT- 
/ (980)^+, / (980) -» vr+vr- 
Non-resonant 



XcO 



K+, 



XcO 



IT ' TT 



0.046 ± 0.030 ± 0.02 



-0.14 ±0.08 ±0.02 



+0.03 
-0.07 



±0.06 ± 0.05 ± 0.02±g;^ 
±0.27 ± 0.12 ± 0.02±g;g§ 
-0.13 ±0.11 ±0.02^;^ 
-0.04 ± 0.08 ± ().i)2-;j; ( ',;| 
-0.30 ±0.22 ±0.02+^;^" 



References 

1. See for example: G. Eilam, M. Gronau and R.R. Mendel, Phys. Rev. Lett. 74, 4984 (1995); 
S. Fajfer, R.J. Oakes, and T.N. Pham, Phys. Lett. B 539, 67 (2002); N.G. Deshpande, 
N. Sinha and R. Sinha, Phys. Rev. Lett. 90, 061802 (2003); M. Gronau, Phys. Rev. Lett. 
91, 139101 (2003); T. Gershon and M. Hazumi, Phys. Lett. B 596, 163 (2004). 

2. A. Abashian et al, Nucl. Instrum. Methods A 479, 117 (2002). 

3. A. Garmash et al. (Belle Collaboration), hep-ex/0412066 



